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ABSTRACT 

Different scaling and autocorrelation characteristics and their application to astro- 
nomical images are discussed: the structure function, the autocorrelation function, 
Fourier spectra and wavelet spectra. The choice of the mathematical tool is of great 
importance for the scaling analysis of images. The structure function, for example, 
cannot resolve scales which are close to the dominating large-scale structures and can 
lead to the wrong interpretation that a continuous range of scales with a power law 
exists. The traditional Fourier technique, applied to real data, gives very spiky spec- 
tra, in which the separation of real maxima and high harmonics can be difficult. We 
recommend as the optimal tool the wavelet spectrum with a suitable choice of the 
analysing wavelet. We introduce the wavelet cross- correlation function which enables 
to study the correlation between images as a function of scale. The cross-correlation 
coefficient strongly depends on the scale. The classical cross-correlation coefficient can 
be misleading if a bright, extended central region or an extended disk exists in the 
galactic images. 

An analysis of the scaling and cross-correlation characteristics of 9 optical and 
radio maps of the nearby spiral galaxy NGC 6946 is presented. The wavelet spectra 
allow to separate structures on different scales like spiral arms and diffuse extended 
emission. Only the images of thermal radio emission and Ha emission give indica- 
tions of 3-dimensional Kolmogorov-type turbulence on the smallest resolved scales 
(160-800 pc). The cross-correlations between the images of NGC 6946 show strong 
similarities between the images of total radio emission, red light and mid-infrared 
dust emission on all scales. The best correlation is found between total radio emission 
and dust emission. Thermal radio continuum and Ha emission are best correlated 
on a scale of about 1' ~ 1.6 kpc, the typical width of a spiral arm. On a similar 
scale, the images of polarised radio and Ha emission are anticorrelated, which remains 
undetected with classical cross-correlation analysis. 

Key words: Physical processes: turbulence - methods: data analysis - galaxies: ISM 
- galaxies: spiral - galaxies: individual: NGC6946 



1 INTRODUCTION 

The increasing resolution of optical, radio, infrared and X- 
ray telescopes over the past decades has greatly improved 
the details visible in astronomical objects in many wave- 
length ranges. The large variety in the structure of external 
galaxies, which have bars, bright central regions, spiral arms 
with clumpy star-forming regions, long known from optical 
photographs, can now also be identified in radio. X-ray and 
IR continuum, and in spectral lines of ionised (Hii), atomic 
(Hi) and molecular (CO) gas components. How these dif- 
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ferent constituents in a galaxy are related is a question to 
which continuous study has been devoted. 

The detection of similar structures in widely separated 
energy ranges gives important information on physical pro- 
cesses and their interplay. For example, the degree of corre- 
lation between the Ha emission and the free-free emission 
in radio continuum tells us about the amount of absorp- 
tion in the optical range. The similarity of far-infrared and 
radio continuum emission components (e.g. in the nearby 
spiral galaxy M 31) indicates that magnetic fields are not 
anchored in the warm medium, but in cool gas clouds (Ho- 
ernes, Berkhuijsen & Xu 1998). Also emission in the CO line 
is correlated with radio continuum emission in spiral arms of 
M31, but it is anticorrelated with the polarised radio emis- 
sion (Berkhuijsen, Bajaja & Beck 1993). 
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Anticorrelations are possibly even more interesting than 
correlations. In the spiral galaxy NGC 6946 a striking anti- 
correlation between the optical spiral arms and the "mag- 
netic arms" seen in radio polarisation was discovered by 
Beck & Hoernes (1996). Although obvious to the eye, the 
analysis of this phenomenon needs sophisticated techniques. 
Based on wavelet analysis, Frick et al. (2000) showed that 
the "magnetic arms" are phase-shifted images of the opti- 
cal arms. Dynamo action is able to generate such structures 
(Rohde, Beck & Elstner 1999). 

The existence of a correlation, anticorrelation or non- 
correlation has important consequences for the interpreta- 
tion of observable quantities which emerge from a combi- 
nation of physical quantities. For example, Faraday rota- 
tion [RM) of polarised radio waves is due to the product 
of electron density (nc) and the magnetic field component 
along the line of sight (-By), averaged along the line of sight. 
Knowledge of < no > (e.g. from dispersion measures of pul- 
sars) is not sufficient to determine the strength of < By > 
unless it is known whether rie and By are correlated, anticor- 
related or non-correlated. If the anticorrelation observed on 
the scale of spiral arms holds also on small scales, the field 
strengths derived from Faraday rotation data of pulsars are 
too small (Beck 2001). 

Relations between images are generally analysed with 
the pixel-to-pixel correlation function. However, this method 
gives little information in the case of an anticorrelation on 
the scale of spiral arms, like in NGC 6946, or when the dif- 
fuse emission on larger scales (e.g. in total radio emission) 
has no counterpart in the other image (e.g. polarised radio 
emission) . 

Another interesting subject is turbulence in the in- 
terstellar medium of galaxies, which has been studied by 
analysing structure functions of appropriate images. Obser- 
vations of radio-wave scintillations in the Milky Way have 
revealed the existence of density irregularities in the diffuse 
ionised gas with scale sizes of > 10^* cm (0.3 pc) down to 
10^ cm. The spectrum of the irregularities is consistent with 
the spectrum of three-dimensional Kolmogorov turbulence 
(Spangler 1999). A spectrum with the same slope was found 
by Minter & Spangler (1996) for both the turbulent fluctu- 
ations in electron density and in Faraday rotation measure 
observed in the same area in the Milky Way. The outer scale 
of these fluctuations is 4 pc. For larger scales up to about 
80 pc Minter & Spangler found a flatter spectrum. As turbu- 
lent scales larger than about ICQ pc are difficult to observe 
in our Galaxy, it is still unknown whether the Kolmogorov 
spectrum extends to large scales. The analysis of images of 
external galaxies may give us this information. 

In laboratory experiments extremely long data series 
are used to get reliable statistics for the slopes of struc- 
ture functions of high order in the inertial range (Monin 
& Yaglom 1971, 1975; Frisch 1995). The analysis of turbu- 
lence in astronomical objects requires observational data of 
very high quality and high angular resolution. The excel- 
lent data set on the external galaxy NGC 6946 seems well 
suited for a study of turbulence on large scales. In a first 
attempt Beck, Berkhuijsen & Uyaniker (1999) analysed the 
high-resolution radio continuum maps of this galaxy using 
the structure function. They found spectra consistent with 
two-dimensional turbulence in the radio continuum emission 
on scales between 0.'5 (twice the beam size) and 5.'0. In this 



paper, however, we show that in the case of images with 
pronounced structures on scales less than 20 x the size of 
the telescope beam, like the bright central area and spiral 
arms covering most of the available images, the structure 
function is not the best analysing tool because it does not 
well resolve these scales. The lack of spectral resolution and 
the relatively poor statistics (i.e. a small ratio of image size 
to beam size) influence the slope of the derived structure 
function, making the interpretation difficult. As astronom- 
ical observations are limited in resolution and extent, the 
mathematical tools used for their analysis have to be chosen 
carefully and the results need a critical evaluation. 

During the last decade a new mathematical tool has 
been developed, wavelet analysis which enables the detec- 
tion of structures of different scales in data sets. First ap- 
plications to astronomical objects showed that wavelets are 
effective in time series analysis (Foster 1996; Frick et al. 
1997a,b), denoising (Tenorio et al. 1999; Chen et al. 2000) 
and structure detection (Frick et al. 2000). 

We have applied two-dimensional wavelet analysis to 
images of NGC 6946 at various wavelengths. Our goal is not 
only to detect the dominant scales in this galaxy, but also to 
see if wavelets can be useful to determine the statistical char- 
acteristics for given scales in the maps. The observation of 
magnetic arms situated in between the gaseous arms, lead- 
ing to an absence of classical cross-correlation between the 
corresponding images in spite of their similarity, stimulated 
this research. 

The paper is organized as follows. In Sect. 2 we present 
the wavelets as a tool of scaling analysis. Their relationship 
with spectra and structure functions is explained in Sect. 3, 
and how wavelets can be used for correlations at a given 
scale is shown in Sect. 4. In the second part of the paper the 
wavelet analysis is applied to images of the nearby galaxy 
NGC 6946 (Sect. 5.1). First its spectral characteristics are 
discussed in Sect. 5.2, then the correlations between pairs of 
images are presented in Sect. 5.3. The results are discussed 
in Sect. 6. 



2 WAVELETS AS A TOOL FOR SCALING 
ANALYSIS 

Wavelet analysis is based on a space-scale decomposition us- 
ing the convolution of the data with a family of self-similar 
basic functions that depend on two parameters, scale and 
location. It can be considered as a generalization of the 
Fourier transformation, which uses harmonic functions as a 
one-parametric functional basis, characterized by frequency, 
or in the case of a space function, by the wavevector k. The 
wavelet transformation also uses oscillatory functions, but 
in contrast to the Fourier transform these functions rapidly 
decay towards infinity. The family of functions is gener- 
ated by dilations and translations of the mother function, 
called the analysing wavelet. This procedure provides self- 
similarity, which distinguishes the wavelet technique from 
the windowed Fourier transformation, where the frequency, 
the width of the window and its position are independent 
parameters. 

We consider the continuous wavelet transform, which in 
the two-dimensional case can be written in the form 
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W(a, x) 



dx 



(1) 



Here x = {x,y), f(x) is a two-dimensional function, for 
which the Fourier transform exists (i.e. square integrated), 
7/) (5) is the analysing wavelet (real or complex, * indicates 
the complex conjugation), a is the scale parameter, and k is 
a normalization parameter which will be discussed below. 

For later considerations the relation between the 
wavelet and the Fourier decomposition will be useful. The 
2-D Fourier transform /(fc) of the function f{x) is defined 



f{x)e '''^dx, 



(2) 



where k = {kx,ky) is the wavevector. Then the inverse 
Fourier transform is 



00 p-\-OC 



}{k)e''"^dk 



00 J —00 



and the wavelet coefficients (pi) can be expressed as 



W{a,x) 



47r2 



f{k)-ip* {ak)e^^ dk. 



(3) 



(4) 



We restrict our analysis to the use of isotropic wavelets. 
It means that the analys ing wave let is an axisymmetric func- 
tion 1/) = P = + y^- The choice of the wavelet 
function depends on the data and on the goals of the anal- 
ysis. For spectral analysis wavelets with good spectral res- 
olution (i.e. well localized in Fourier space, or having many 
oscillations) are preferable, for local structure recognition a 
function, well localized in the physical space, is preferable. 
(Note that the spectral resolution Afc and the space res- 
olution Ax are strongly related and are restricted by the 
uncertainty relation AxAA; > 27r.) An obligatory property 
of the wavelet is the zero mean value J J ip{x,y)dxdy — 0. 

In the case under consideration the choice of the wavelet 
was determined by the wish to have more independent points 
for further structure analysis which led to a simple real 
isotropic wavelet with a minimal number of oscillations, 
known as the Mexican Hat 



(5) 



We shall refer to this function as MH. 

For a better separation of scales (to analyse spectra and 
to find the scale of dominant structures) another isotropic 
wavelet was used which is defined in Fourier space by the 
formula 



TT < |fc| < 47r 
\k\ < TT, \k\ > Atv 



(6) 



The function is localized in Fourier space in a ring with a 
median radius 2tt and vanishes for \k\ < tt and |A;| > 47r. 
This wavelet definition provides a relatively good spectral 
resolution, it will be referred to as PH (it is our Pet Hat). Q 
In physical space the PH wavelet is obtained by numerical 
integration of Both wavelets MH and PH are shown in 
Fig. I 

t This function was introduced in Aurell et al. (1994) for model- 
ing the two-dimensional turbulence. 




Figure 1. Two isotropic wavelet functions used in this paper: the 
MH wavelet (thin line) and the PH wavelet (thick line). 



The wavelet transform (|l|) is unique and reversible 
which means that the analysed function f{x, y) can be re- 
constructed from its wavelet decomposition. In our analysis 
we do not need the inverse transform and do not give the 
reconstruction formula. (An extended description of contin- 
uum wavelet transform can be found in e.g. Holschneider 
(1995) and Torresani (1995).) 

To illustrate how wavelets decompose the image in dif- 
ferent scales we show in Fig. ^ the optical broadband emis- 
sion image of the galaxy NGC 6946 and its wavelet coeffi- 
cients W{a,x) for three different scales a. For this example 
the PH wavelet was used. 



3 SPECTRA AND STRUCTURE FUNCTIONS 

In the studies of scaling properties of turbulence a commonly 
used characteristic of a turbulent field is the spectral energy 
density E{k) which includes the energy F{k) — |/(fc)p of all 
Fourier harmonics with wavenumbers k, for which \k\ — k 

E{k) = / F{k)dk . (7) 

J\k\=k 

The spectral energy is related to the autocorrelation function 
/ fix)f{l- x)dx (8) 

-oo ^ — oo 

by its Fourier transform 

F{k) = C{k) = / / C{l)e-''^'dl . (9) 



oo ^ — oo 



The vector I defines the shift of the image in the convolu- 
tion ^ and is the scale parameter. In the important case 
of isotropic turbulence the autocorrelation function depends 
only on the distance between two points C(Z) = C{1), and 
the spectral energy depends only on the modulus of the 
wavevector F{k) = F{k). Then these two functions are re- 
lated by the Hankel transform 



(10) 



F{k) = 2-K I C{l)Jo{kl)ldl , 



where Jo is the Bessel function and E{k) — 2nkF{k). 

Another often used characteristic for scaling studies is 
the structure function defined for arbitrary order q as 
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Figure 2. NGC6946. Map of the optical emission in red light (central region subtracted) and its wavelet decompositions for 3 different 
scales: 0.'5, 1.'4 and 4'. The size of each image is 11' X 11'. 



s.(o = <(/(s)-/(^-or>|,-|=, , (11) 

where the brackets < ... > mean the average value. Cal- 
culation of high-order structure functions requires a high 
accuracy of the initial data. In the case of maps of external 
galaxies, where a relatively small number of grid points is 
available and the noise is significant, only the second-order 
function &, corresponding to the energy spectrum (^), can 
be discussed. 

In the wavelet representation the scale distribution of 
the energy can be characterized by the wavelet spectrum, 
defined as the energy of the wavelet coefficients of scale a of 
the whole physical plane 

M(a) = / / \W{a,x)\'dx . (12) 



The wavelet spectrum can be related to the Fourier spec- 
trum. Using (^ one can easily rewrite (^2|) in the form 

M(a) = ^ / / \f{kt\^{ak)\'dk . (13) 

This relation shows that the wavelet spectrum is a smoothed 
version of the Fourier spectrum. In the isotropic case ( |l3| ) 
has a more simple form 

M{a) = ^ / E{k)\^lj{ak)fdk . (14) 

8^ Jo 

A generic property of fully developed small-scale turbu- 
lence is the presence of an inertial range of scales in which 
the energetic characteristics E{k), 82(1) and M(o) follow 
power laws. Let us consider the relation between the spec- 
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tral indices of these characteristics. Let the structure func- 
tion follow a power law 

SiiO'^l''. (15) 

Then, F{k) ~ k~'^~^ and the spectral energy density ex- 
hibits a power law 

E{k) ~ fe-(^+'> . (16) 

The behaviour of the wavelet spectrum M(a) depends on the 
normalization k in definition (^. We will use k = 2, which 
gives the same power law for the wavelet spectrum as for 
the structure function 



M{a) 



(17) 



Choosing k = 3/2 the wavelet spectrum becomes M{a) ~ 
a^"*"^, which is convenient if the result should be compared 
with the Fourier spectrum E{k). It should be noted that 
the scale parameter a commonly used in wavelets has the 
same meaning as the distance I in the autocorrelation (^ 
or structure function (pT[). Below we shall normally use the 
character a for the scale parameter in any spectral charac- 
teristic. 

An important remark concerns the power-law search in 
spectra. Every wavelet has its own spectral portrait and the 
tail of the corresponding spectral energy distribution often 
itself follows a power law. The power index of this law defines 
the absolute limit of the spectral slope that can be obtained 
using the given wavelet. 

Note that the traditional calculation of the structure 
function ^2 following (^ij) can be interpreted as the calcula- 
tion of the wavelet spectrum (|l^), obtained using a "special" 
anisotropic wavelet, which gives the difi'erence between two 
delta functions separated by the unit distance 

'ip{x,e) = 5{x)-5{x-e) , (18) 

where e is the unit vector. This quasi-wavelet (|l^) is a bad 
wavelet in the sense that being very well localized in physical 
space it has inevitably a very bad spectral resolution. This 
means that structure functions give a poor scale separation 
and as a result they are always very smooth. This is why 
the structure functions can be useful only in the case of 
a well developed range of scales for which a power law is 
established. 

We conclude this section with an illustrative example 
which is relevant to further data analysis. In Fig. |^ two 
images are presented, the second one being a phase-shifted 
copy of the first one. The spectral properties of these two 
images are identical and are shown in Fig. The upper 
curve shows the structure function S2 (pT|), which is almost 
flat over a large range of scales. A steep decrease appears 
at large scales when the distance I surpasses the total size 
of the structure. The next two curves reproduce the wavelet 
spectra obtained with different wavelets, MH (^) and PH 
(^). The graphs clearly illustrate the increase of spectral 
resolution: some maxima are indicated in the MH spectrum 
and three maxima are well pronounced in the PH spectrum. 
These maxima correspond to the scale of the width of an 
individual spoke in the image, the scale of its length and the 
global scale of the structure. The last (lower) curve shows the 
Fourier spectrum, calculated via the Hankel transform ( [lo| ) 
of the autocorrelation function (^). It displays of course the 
best spectral resolution but includes the upper harmonics - 



r„(a) 




Figure 3. Spectral and correlation characteristics of artificial im- 
ages, a Test images with sides of length 2; b spectral character- 
istics of a single image: second-order structure function (crosses), 
MH wavelet spectrum (black dots), PH wavelet spectrum (cir- 
cles), Fourier spectrum (line) (the spectra are shifted along the 
vertical axis for better presentation), a = 27r/fc; c wavelet cross- 
correlation coefficient. 



two maxima at a < 0.25 are the double harmonics of the 
central structure. Note that the Fourier spectrum is shown 
versus the physical scale a which can be interpreted as the 
wavelength, a — 2n/k, thus the high frequencies are to the 
left (small a). 

This example illustrates the usual problem of Fourier 
presentation of observational data - real data give very spiky 
spectra and the interpretation of different maxima in them is 
difficult. A reliable choice of wavelets leads to smooth spec- 
tra which nevertheless clearly indicate the dominant scales 
in the analysed data. Note also that the global scale indi- 
cated by the structure function and the wavelets at a « 3 is 
not visible in the Fourier spectrum - the size of the image 
is too small to contain the corresponding Fourier harmonic. 
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4 WAVELET CROSS-CORRELATIONS 

In this section we discuss correlations of images of the same 
extragalactic object at different wavelengths and how they 
can be quantified. 

Let us consider two images (maps), y) and f2{x, y), 
with the same angular resolution and the same number of 
pixels. The simplest way for a linear correlation study is the 
direct calculation of the correlation coefficient point by point 
(pixel by pixel) 



< /i >)(/2»- < h >) 



(19) 



The accuracy of this estimate depends on the degree of cor- 
relation and on the number of independent points n (see for 
example Edwards (1979)) 



Ar = 



(20) 



Here n is the ratio of the map area to the area of the beam 
(n ~ 1700 in the maps discussed in Section 5). 

The "pixel by pixel" correlation coefficient is a global 
characteristic which contains all the scales present in the 
images, including the largest one which reflects the fact that 
the centre of the extended intensity distribution is nearly at 
the same position at most wavelengths (see the last panel of 
Fig. i). 

It is interesting to investigate whether the correlation 
coefficient depends on scale. Therefore we introduce the cor- 
relation coefficient for a given scale a 



,(a) = 



/ / W^_{a,x)Wi(a,x)dx 



(21) 



(Afi(a)M2(a))i/2 

To estimate the error Ar^, (a) we use formula ( po| ) taking 
n = {L/af , where L is the size of the map. 

To illustrate the behaviour of this wavelet correlation 
(i.e. correlation on agiven scale) ^ let us return to the ex- 
ample given in Fig. |3[ The two images shown in the upper 
panel have an angular shift tt/IO (the spokes in the second 
map are just in between those in the first one). Calculating 
the "pixel by pixel" correlation coefficient yields rp — —0.01. 
This vanishing value of the correlation coefficient reflects the 
fact that the bright details of both images occupy a small 
part of the image and do not overlap. At the same time it 
is clear that the images are very well correlated in general, 
because the images present two identical structures, one of 
which is rotated. In Fig. ^ we present the wavelet correla- 
tion coefficient ru,(a) calculated for these two images. The 
plot explicitly shows that only the smallest scales are re- 
ally non-correlated. On the scale corresponding to the mean 
distance between the spokes a strong anticorrelation is vis- 
ible (r™ — —0.75). Note that the minimum in r^, does not 
coincide with the maximum in the spectra - the latter corre- 
sponds to the width of an individual spoke, and the former to 
the distance between them. On the scales which characterize 
the length of the spokes, the images become correlated and 
tends to unity. On larger scales the pictures are identical. 

The relationship between the correlation coefficient rp 



t The wavelet correlation was first introduced by Nesme-Ribes 
et al. (1995) for two time series of solar data. 



and the wavelet correlation can be obtained using (^, 
(p^, (^), and in the case under consideration (k = 2) has 
the form 

/ r^{a){Mi{a)M2{a)y^^a-'^da 
^ {J Mi{a)a-^da J M2(a)a-ida)i/2 ' 

We conclude this section by a comparison of the intro- 
duced wavelet correlation with the cross-correlation func- 
tion, defined as 



(22) 



Hi)= / / Mx)f2{l-x)dx . (23) 

Using (§),(§|),(|l^) and (^l|) one gets the following relation 
/ / h(k)\-4,(ak)\'^dk 



rw{a) 



(/ / Ci{k)\iP{ak)\'2dk J J C'2(fc)|7/i(afc)|2dfc)i/2 



(24) 



which in the limit of a wavelet with ideal scale resolution 
(wavelet tends to become a harmonic) leads to a relation 
which in the isotropic case can be written as 



,{a) = 



h{2Tv/a) 



(Ci(27r/a)C2(27r/a))i/2 



(25) 



We emphasize that all the relations between wavelet co- 
efficients, Fourier decomposition and correlation functions 
are exact only if the limits of the integrals surpass the real 
range of scales (spatial frequencies) present in the object. 
From the practical point of view, in the case of noisy data 
with relatively poor spatial resolution, the calculation of 
ru,(a) via wavelets has a strong advantage because in this 
way we first separate the structures of a given scale present 
in the data. The correct calculation of the Fourier trans- 
form h{2n/a), however, implies accurate knowledge of the 
cross-correlation function (^3|) in the whole range of scales 
(the same is true for the spectral energy density and the 
autocorrelation function). 

In the simple case of the test in Fig. ^, both wavelets 
(MH and PH) give practically the same correlation function 
rw{a). In this purely symmetric example the same result 
can be obtained using Eq. (psj). The Fourier decomposition 
was used by Elmegreen et al. (1992) to search for simple 
symmetries in galactic images. However, the symmetry in 
such maps is generally weak and superimposed onto a noisy 
background. Reliable statistics will then be important to get 
significant results. Therefore we prefer the wavelet with the 
best spatial localization and we use the MH wavelet for the 
correlation study which occupies a smaller area than the PH 
wavelet for the same scale. We note again that every Fourier 
harmonic covers the whole plane so that its spatial location 
is completely absent. 



5 ANALYSIS OF SCALES AND STRUCTURES 
IN OPTICAL, RADIO AND INFRARED 
MAPS OF NGC 6946 

5.1 The data 

NGC 6946 is a nearby spiral galaxy for which a distance of 
5.5 Mpc is usually assumed. Data in many spectral ranges 
are available. For our analysis we used: 

a) The map of the optical broadband emission in red 
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Figure 4. Images of NGC6946 in different spectral ranges. Top row (from left to right): a optical (red) light, b linearly polarised radio 
continuum emission at A6.2 cm, c total radio continuum emission at A6.2 cm. Middle row: d nonthermal radio continuum emission at 
A3. 5 cm, e thermal radio continuum emission at A3. 5 cm, f Ha line emission of ionised hydrogen Hll. Bottom row: g line emission 
of neutral atomic hydrogen Hi, h (1-0) line emission of the CO molecule, i mid-infrared emission of warm dust and PAH particles at 
12—18 /im. The size of every image is 11' X 11' and the beam diameter is always 15". 



light (R) from the digitized Palomar Sky Survey.^ We sub- 

§ Based on photographic data of the National Geographic Society 
— Palomar Observatory Sky Survey (NGS— POSS) obtained using 
the Oschin Telescope on Palomar Mountain. The NGS— POSS was 
funded by a grant from the National Geographic Society to the 
California Institute of Technology. The plates were processed into 
the present compressed digital form with their permission. The 
Digitized Sky Survey was produced at the Space Telescope Sci- 
ence Institute under US Government grant NAG W— 2166. 



traded foreground stars from the map and smoothed it to 
the angular resolution of the radio continuum maps. The 
result is shown in Fig. ^. 

b) The map of the linearly polarised radio continuum 
emission at A6.2 cm (PI6), combined from observations with 

the VLA synthesis telescope operated by the NRAO|^ and 
the Effelsberg single dish operated by the MPIfR (Beck & 

^ The National Radio Astronomy Observatory is a facility of the 
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Hoernes 1996). The angular resolution is 15" which cor- 
responds to 400 pc at the assumed distance to NGC 6946 
(5.5 Mpc). 

c) The map ol the total radio continuum emission at 
A6.2 cm (16), combined from observations with the VLA 
synthesis telescope and the EfFelsberg single dish (Beck & 
Hoernes 1996). The angular resolution is 15". We subtracted 
a strong unresolved source on the northwestern edge of 
NGC 6946 because it seems unrelated to the galaxy. 

d) The map of the nonthermal (synchrotron) radio con- 
tinuum emission at A3. 5 cm (NTH3), derived from the spec- 
tral index map between A3. 5 cm and A20.5 cm at 15" res- 
olution (Beck, in prep.) assuming a constant synchrotron 
spectral index of a = —0.9 (where the flux density S varies 
with frequency v as Sv oi 1/°'). The spectral index and non- 
thermal intensity were computed only for pixels where both 
intensity values exceed 10 times the rms noise. 

e) The map of the thermal radio continuum emission 
(free-free emission from ionised hydrogen) at A3. 5 cm (TH3) 
which is the difference map between the total and the syn- 
chrotron emission. 

f) The map of the Ha line emission (Ha) of ionised 
hydrogen H 11, integrated over the whole frequency width of 
the line (Ferguson et al. 1998). We subtracted foreground 
stars and smoothed the map to the angular resolution of 
the radio continuum maps. The Ha map by Ferguson et al. 
(1998) has a higher signal-to-noise ratio than the map used 
in Frick et al. (2000). 

g) The map of the radio line emission of neutral atomic 
hydrogen Hi at A21.1 cm (Hi) observed with the Wester- 
bork Synthesis Radio Telescope, integrated over the whole 
frequency width of the line (Kamphuis & Sancisi 1993). The 
Hi map has an original resolution of 13" x 16" and was not 
smoothed. 

h) The map of the (1-0) line emission of the molecule 
CO at A2.6 mm (CO) observed with the IRAM 30-m single- 
dish telescope, integrated over the whole frequency width of 
the line (Walsh et al. 2001). The CO map has an original res- 
olution of 22" and was not smoothed. The lower resolution 
of the CO map compared with the other maps may influence 
the comparison on small scales, but the effect turns out to 
be small (see Fig. p^ ). 

i) The map of the mid- infrared emission (12-18 ^m) 
(DUST) observed with the ISOCAM camera (filter LW3) on 
board of the ISO satellite (Dale et al. 2000) . This wavelength 
range is dominated by the band of emission lines from poly- 
cyclic aromatic hydrocarbons (PAHs) at 12.7 /xm, but also 
includes continuum emission from small, warm dust grains. 
The LW3 map has an original resolution of about 8" and 
was smoothed to a beamwidth of 15". 

The map in LW2 (5-8.5 jim) is fully dominated by PAH 
emission bands. As it is very similar to the LW3 map, we 
did not include it in our analysis. 

All maps were transformed to the same area (11' x 11') 
and pixel size (2"). Each map contains 331 x 331 pixels. Plots 
of these images are shown in Fig. ^. 

In the radio continuum maps, we removed the bright 
double source at the western edge of the galaxy which has 
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no counterpart in any other spectral range and thus is most 
probably unrelated to NGC 6946. Finally, we removed the 
central region in 7 of the 9 maps by subtracting fitted Gaus- 
sian sources with half-power widths of about 0.'3 (infrared, 
total, nonthermal and thermal radio), 0'5 (Ha, CO) and 1' 
(red light). This avoids artefacts produced by improper con- 
tinuum background subtraction in the central region of the 
Ha map as well as contamination of the mid-infrared map 
by enhanced continuum emission from warm dust grains. 

5.2 Spectral characteristics 

Before starting the presentation of spectral characteristics 
of maps of an external galaxy we want to emphasize that 
one should be very careful when interpreting such kind of 
results for two reasons. Firstly, for the scales that are reli- 
ably resolved, the maps cannot be treated as homogeneous 
and isotropically turbulent (as is possible when analysing 
the Milky Way turbulence from Galactic surveys (Minter & 
Spangler 1996)). One can expect to find the dominant scale 
in the map at a certain wavelength (if it exists) and the scale 
below which the field under consideration can be interpreted 
as homogeneous - that is one may separate the large-scale 
structures and the small-scale turbulence. Secondly, one has 
to carefully separate the spectral properties of the analysed 
image from the spectral properties of the analysing func- 
tion and from structures caused by observational problems 
(regions without observations, bright sources, instrumental 
noise etc). 

We first show the spectral characteristics as defined 
above for two maps of NGC 6946: the total radio emission 
(Fig. 1^ and the distribution of the Ha emission (Fig. 
Clearly, the structure function totally smoothes the spectral 
energy distribution suggesting a spectrum of fully developed 
turbulence on all scales with a power law of slope 0.7. How- 
ever, the other curves show that such a spectrum is not 
present in the data. The Fourier spectrum clearly displays 
the border between the informative part of the spectrum and 
the small-scale noise, visible as a flat (or even increasing, as 
in the Ha spectrum) part for scales smaller than 0.'2. This 
scale is close to the size of the beam (0^25) and proves the 
evident statement that little can be said about scales below 
the scale of the beamwidth. The structure function and the 
wavelet spectra show less flattening on scales smaller than 
the beam. 

Wavelet spectra again give an intermediate presentation 
of the spectral properties. As the optimal spectral charac- 
teristic the PH wavelet spectrum can be proposed which 
smoothes the numerous peaks in the Fourier decomposition 
but is still sensitive enough to indicate the scales contain- 
ing most of the energy. The 16 spectrum shows two weak 
maxima at 0^5 and l!3. Note that the spectrum obtained 
by the MH wavelet does not separate these two maxima. 
The wavelet spectrum of Ha shows a well pronounced maxi- 
mum at 0!8-lf0. Even this maximum is hardly visible in the 
structure function. 

In Fig. ^ we present wavelet spectra calculated for all 9 
images of NGC 6946 using the PH wavelet (^. The range of 
scales is limited to 0.'2 < a < 4^0. Comparing the spectral 
distributions we note that only Ha and TH3 show a clear 
maximum (though shifted in scale, see Sect. 6.1). Less strong 
maxima are visible in CO, Hi, R and nonthermal emission 
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Figure 5. NGC6946. Spectral characteristics of the map of total 
radio emission (16) ; second-order structure function (crosses) , MH 
wavelet spectrum (dots), PH wavelet spectrum (circles), Fourier 
spectrum (line). The size of a pixel is 0.'03 and the beam diameter 
is 0.'25. 
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Figure 6. NGC 6946. Spectral characteristics of the Ha map: 
second-order structure function (crosses), MH wavelet spectrum 
(dots), PH wavelet spectrum (circles), Fourier spectrum (line). 
The size of a pixel is 0.'03 and the beam diameter is 0.'25. 



NTH3. Almost flat are the spectra of total radio emission 
16 and polarised radio emission PI6 - they both show weak 
maxima near O'.b and l'.2. The location of the maxima differs 
for the images. The thermal emission TH3 and Hi contain 
most of the energy on the scale of about O'.b, CO on 0'.7, 
Ha and R on l!0, the NTH3 (and maybe 16) on l!4-l.'5. 
For comparison: the typical size of a giant complex of gas 
clouds is « 1 kpc (0.'6), the typical width of the spiral arms 
is X 1.5 kpc (0.'9). 

The spectra in Fig. |^ show that the spatial resolution of 
the available images is insufficient for the study of properties 
of the small-scale turbulence. The inertial range is expected 
to start at scales several times smaller than the dominant 
scale (turbulent macroscale). If the observed object contains 
the classical three-dimensional Kolmogorov turbulence with 
the spectrum E{k) ~ k"^^'' (corresponding to the structure 
function S2 ~ l^^^), the two-dimensional projection should 
be characterized be a structure function S2 ~ l^^^ (Spangler 
1991). The spectral slope of the projected field steepens be- 
cause the smaller scales are smoothed in the projection. In 
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Figure 7. NGC 6946. Wavelet spectra (analysing wavelet is PH). 
The thin lines above Ho and TH3 indicate a slope of 5/3. 
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Figure 8. NGC 6946. Scale-by-scale wavelet cross-correlations: 
Ha-PI6 (dots); Ha-TH3 (stars); Ha-NTH3 (circles). 

Fig. ^, only TH3 and Ha show Kolmogorov-type behaviour 
on small scales (see Sect. 6.1). 

5.3 Cross-correlations in the galaxy NGC 6946 

We start by calculating the classical "pixel by pixel" cross- 
correlation coefficient rp for every pair of galactic maps un- 
der discussion. The values of rp are given in the upper part 
(above the diagonal) of Table 1. Note that for most pairs 
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the value of the coefficient indicates a statistically signifi- 
cant correlation. Only the PI6 map in pair with TH3, Ha 
and H i yields correlation coefficients significantly less than 
0.5 (0.23, 0.25 and 0.29, respectively). 

For each pair of maps the wavelet correlation coeffi- 
cient was also calculated, using the MH wavelet. As an 
example, the behaviour of (a) is shown in Fig. ^ for three 
pairs of data: Ha-PI6, Ha-TH3 and Ha-NTH3. The high 
correlation on largest scales reflects only the coincidence of 
the extended regions corresponding to the galaxy's disk (see 
Fig. 1^, bottom right). However, on intermediate scales the 
three curves display very different behaviours. The anti- 
correlation between Hq and PI6 is of special interest (see 
Sect. 6.2). The tendency of all curves in Fig. | to grow to- 
wards the smallest scale (a = 0^125) indicates that this scale 
is affected by fluctuations on scales significantly smaller than 
the resolution of the maps (0'25) and thus will not be con- 
sidered further, whereas the next computed scale {a — 0'.2), 
though still slightly smaller than the beam size, gives useful 
information. Note that the number of independent points 
decreases towards larger scales which results in increasing 
error bars on these scales. 

To illustrate the difference in correlating the initial data 
and the wavelet decomposition on various scales, we show 
the corresponding plots in Fig. ^ for two extreme cases of 
Fig. ^ The correlation of the maps PI6 and Hq gives rp = 
0.25. This low value of rp is understandable when looking at 
the distribution of points in Fig. ^a. In panel (b) we show the 
wavelet coefficients W{a, x) obtained for the PI6 map versus 
the corresponding wavelet coefficients calculated for the Ha 
map. The scale parameter is fixed at the value a = l!0, which 
corresponds to the minimum in the curve (a) for this pair 
of maps in Fig. |^. The cloud of points becomes elongated 
along the line y = —x which corresponds to a negative value 
of r^,. Note that the wavelet coefficients W{a,x) become 
negative in case of local minima on a certain scale in the 
analysed data. Panel (c) shows the corresponding plot for 
the same pair of maps, but for the scale a — 2.'8, which gives 
rw = 0.57. The tendency of a positive correlation becomes 
visible. 

Panels (d), (e) and (f) of Fig. ^present the analogous 
plots for the maps TH3 and Hq. This pair of data shows 
a high degree of correlation: rp = 0.81. Note that the plot 
in panel (e) corresponding to W{a,x) at a = l.'O (a local 
maximum in r^ia), see Fig. ^) gives a better confined cloud 
of points than the initial data in panel (d), in spite of the 
fact that the correlation coefficient is r^ = 0.71 (less than 
rp). This means that the scatter in the points in panel (d) 
is compensated by the larger number of points along the 
central strip of the cloud. The points in panel (e) are also 
more homogeneously distributed than in panel (d). Panel 
(f) shows again the scale a — 2.'8, on which TH3 and Hq 
are less correlated than on the scale a — l.'O (~ 1.6 kpc, the 
typical width of a spiral arm). 

When analysing galactic images, the scales which char- 
acterize the width of the spiral arms are of special interest. 
Therefore we calculated the mean value of r^ for the range 
of scales 0.'4 < a < l'.2. These scales are also more certain in 
the sense that they are large enough with respect to the size 
of the beam and small enough to give satisfactory statistics 
over the area of the galactic image. The values derived are 
given in the lower part of the Table 1 (below the diagonal) . 



Comparing the corresponding values of rp and r„ in 
Table 1, three groups of map pairs can be distinguished: 
1) the two correlation coefficients are similar (the difference 
between rp and r^ is less than 25%), 2) the two correlation 
coefficients strongly differ, rp/r^ ~ 2 or even r,„ < {rp is 
never negative in the table), 3) comme ci - comme ga. 

Group 1 includes most of the pairs containing Hi or 
TH3 (only the pairs H1-PI6, TH3-PI6 strongly drop out 
of this rule) and several pairs with Hq: Hq-R, Hq-I, Hq- 
DUST. The coefficients rp and practically equal for 

H1-I6, Hi-NTHS, TH3-HQ, and TH3-I6. In the whole table 
only the pairs Hi-CO and Hi-DUST give r„ > rp. 

Group 2 includes first of all the correlations with PI6. 
Three pairs (PI6-R, PI6-TH3, PI6--Hq) yield negative r„. 
Strong differences between r^ and rp occur in R-NTH3 
(0.29/0.86), R-CO (0.35/0.83), NTH3-C0 (0.32/0.80), 
NTHS-Hq (0.34/0.64), and NTH3-DUST (0.47/0.89). 

We calculated 36 cross-correlation functions r^i^a), one 
for each pair of galactic maps. Because of space limitation 
we present the correlation curves r^i^a) only for two sets of 
pairs (see Fig. ^ red emission R (left column) and polarised 
emission PI6 (right column) correlated with each of the other 
8 maps. 

For most pairs r^ varies considerably with scale. Only 
a few curves display a weak dependence on scale, for ex- 
ample R-DUST, which leads to similar values for r^ and 
rp (Group 1). The best coincidence of r„ and rp, however, 
is obtained when r„ has a maximum on the scales chosen 
for the calculation of r^ for Table 1. The pairs R-Hi (see 
Fig. |l^) and Hq-TH3 (Fig. |) are good examples of this 
case. Also the pairs Hi-CO, Hq-I6, TH3-DUST show weU 
pronounced local maxima on the scale of the arm width. 

Obviously, the main part of Group 2 (high contrast be- 
tween rp and r^) consists of pairs for which the r„ curve 
has a minimum on intermediate scales. Apart from all pairs 
including PI6, the pairs CO-NTH and CO-R belong to this 
group. In these cases the classical cross-correlation coeffi- 
cient can be misleading (see Sect. 6.2). 



6 CONCLUSIONS AND DISCUSSION 

In this paper we have used the wavelet analysis to find the 
scales of dominant structures in maps of the spiral galaxy 
NGC 6946 observed in various spectral ranges and to study 
cross-correlation coefficients between pairs of maps. Our 
conclusions are summarised below. 



6.1 Wavelet spectra 

We would like to stress once again that the scaling analy- 
sis of observational data requires an adequate choice of the 
mathematical tool. We have shown that the structure func- 
tion, which is a commonly used method for scaling analysis 
of small-scale turbulence, is not suitable for scales which are 
close to the dominating large-scale structures (see Figs. ^ 
and Due to poor scale resolution, the structure function 
may suggest that a continuous range of scales with a power 
law exists (as e.g. in Beck et al. 1999) even when this is not 
present in the data. On the contrary, the traditional Fourier 
technique applied to real data gives very spiky spectra in 
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Table 1. Cross-correlations between 9 maps of NGC6946. The classical correlation coefficients Vp are given above the diagonal of the 
table. Below the diagonal the mean values of the wavelet correlation coefficient ru,(a), calculated for the range of scales 0f4 < a < 1.'2, 
are shown. 



R PI6 16 NTH3 TH3 Ha Hi CO DUST 



R 




0.51 
±0.03 


0.84 
±0.02 


0.86 
±0.02 


0.66 
±0.03 


0.77 
±0.02 


0.44 
±0.03 


0.83 
±0.02 


0.88 
±0.02 


PI6 


-0.11 
±0.06 




0.56 
±0.03 


0.67 
±0.03 


0.23 
±0.03 


0.25 
±0.03 


0.29 
±0.03 


0.53 
±0.03 


0.49 
±0.03 


16 


0.49 
±0.05 


0.19 
±0.06 




0.92 
±0.02 


0.81 
±0.02 


0.77 
±0.02 


0.47 
±0.03 


0.78 
±0.02 


0.94 
±0.02 


NTH3 


0.29 
±0.06 


0.36 
±0.06 


0.71 
±0.04 




0.60 
±0.03 


0.64 
±0.03 


0.40 
±0.03 


0.80 
±0.02 


0.89 
±0.02 


TH3 


0.57 
±0.05 


-0.10 
±0.06 


0.77 
±0.04 


0.53 
±0.05 




0.81 
±0.02 


0.45 
±0.03 


0.57 
±0.03 


0.79 
±0.02 


Ha 


0.68 
±0.05 


-0.22 
±0.06 


0.65 
±0.05 


0.34 
±0.06 


0.80 
±0.03 




0.50 
±0.03 


0.65 
±0.03 


0.81 
±0.02 


Hi 


0.35 
±0.06 


0.14 
±0.06 


0.47 
±0.06 


0.39 
±0.06 


0.37 
±0.06 


0.41 
±0.06 




0.41 
±0.03 


0.43 
±0.03 


CO 


0.35 
±0.06 


0.24 
±0.06 


0.48 
±0.06 


0.32 
±0.06 


0.36 
±0.06 


0.36 
±0.06 


0.49 
±0.06 




0.85 
±0.02 


DUST 


0.56 
±0.05 


0.14 
±0.06 


0.84 
±0.04 


0.47 
±0.06 


0.63 
±0.05 


0.60 
±0.05 


0.57 
±0.05 


0.66 
±0.05 





which the separation of real maxima and high harmonics 
can be difficult. 

We have shown that wavelets (provided the choice of 
the wavelet is adequate to the aim of the analysis and to the 
analysed data) allow to obtain reliable spectral characteris- 
tics of the object even if the range of scales in the map is 
relatively small (about one and a half decades in our cases) . 
A suitable choice of wavelets leads to smooth spectra which 
clearly indicate the dominant scales in the analysed data. 

Using wavelets, the spectral characteristics of 9 optical, 
radio and infrared maps of the galaxy NGC 6946 were cal- 
culated. The wavelet spectra in Fig. |^ show three types of 
curves: 

a) R, 16, CO and DUST increase smoothly with increas- 
ing scale. These maps show structures on small and large 
scales (clouds, spiral arms, an extended central region and 
extended diffuse emission). 

b) The spectrum of PI6 is practically flat on all scales 
because this map neither contains an extended central region 
nor extended diffuse emission. 

c) TH3 and Ha show a prominent maximum on inter- 
mediate scales corresponding to the width of the spiral arms. 
The maximum in TH3 occurs on a somewhat smaller scale 
than that in Ha. A possible explanation is the increase of 
absorption of the Hq emission on small scales where gas 
densities are generally higher (see also Fig. [l^ ). 

We conclude that the wavelet spectrum reveals signa- 
tures of different phenomena in an object and allows to com- 
pare their relative importance within the object as well as 
between different objects or wavelength ranges. 

Only two spectra, Hq and thermal emission TH3, in- 



clude an increasing part on small scales, which (with great 
caution) can be interpreted as an indication of the Kol- 
mogorov inertial range in which the energy is transferred 
from large to small scales. In Fig. ^ we show the expected 
theoretical slope of 5/3 as thin lines at a < 0.'5 (800 pc) 
for Ha and at a < 0.'35 (560 pc) for TH3. We also anal- 
ysed the spectrum of the Ha image at full resolution (~ 2") 
and found that the slope of 5/3 continues down to a scale 
of about O'l (160 pc). Note that this slope should charac- 
terize the two-dimensional projection of a three-dimensional 
turbulent field with Kolmogorov spectral properties. 

Some spectra in Fig. Q include an increasing part with 
a slope of less than 1. Red light R shows a slope ~ 0.6 for 
0:25 < a < 1', in CO the slope is close to 0.7 for 0^25 < 
a < 1', and in the same range of scales the DUST emission 
is characterized by a slope « 0.5. The total radio emission 
16 shows a large range of scales with growing energy, but 
the slope is not constant and depends on the choice of the 
range; a weak maximum near a — O'.b may indicate two 
subintervals with different slopes. The mean slope in the 
range 0^25 < a < I'A is about 0.6-0.7. 

Structure functions with a slope of about 2/3 over a 
range of scales have been interpreted as an indication of two- 
dimensional turbulence (Minter & Spangler 1996; Beck et al. 
1999). This interpretation, however, is uncertain because the 
physical behaviour of two-dimensional turbulence is quite 
different from three-dimensional turbulence. 

In 2-D turbulence two inertial ranges are possible 
(Kraichnan 1967). This is due to an additional conserva- 
tion law, i.e. in the limit of small viscosity the total square 
of vorticity, called enstrophy, becomes an integral of mo- 
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Figure 9. NGC6946. Cross-correlations of the intensities of the initial images and of their wavelet decompositions. Polarised radio 
emission (PIG) versus Ha: a initial data, b wavelet coefficients w for scales a = 1', c for a = 2f8. Thermal radio emission {TH3) versus 
Ha: d initial data, e wavelet coefficients w for scales a = 1', f for a = 2.'8. 



tion as well as the total energy. The first inertial range, the 
range of inertial transport of enstrophy to smaller scales, 
appears on scales smaller than the scale where the turbu- 
lence is excited. The dimensional arguments lead to a very 
steep spectrum E{k) ~ , which corresponds to a struc- 
ture function S2 ~ ■ No average spectral flux of energy 
is expected in this inertial range. The second inertial range 
can appear on scales larger than the scales on which the 
energy input into the flow occurs. Although characterized 
by the usual Kolmogorov spectral slope {E(k) ~ k~'^^^ or 
S2 ~ l^'^), it has very different physical properties because 
the energy in this inertial range can be transferred only to 



larger scales (so-called inverse energy cascade). It is impor- 
tant to note that the formation of the inverse energy cascade 
is a relatively slow process and implies a stationary forcing 
on intermediate scales (see, e.g., Borue (1994) or Babiano, 
DubruUe & Frick (1997)). Two-dimensional turbulence with 
a slope of 2/3 implies the existence of an inertial range of 
scales with inverse cascade. This means that the energy of 
turbulent motion on these scales does not emerge from the 
larger (galactic) scale, but is transferred from smaller (sub- 
galactic) scales on which nevertheless the motion already 
has to be two-dimensional, which is improbable. 

At the moment we cannot explain the origin of the 
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Figure 10. NGC6946. Wavelet cross-correlations for red light (R, left) and polarised radio emission (PI6, right) versus each of the other 
maps. 



slopes of about 2/3 detected in the wavelet spectrum of R, 
CO, DUST and 16 on small scales. Furthermore, we do not 
understand why a slope of about 5/3 is observed only for the 
ionised gas (Ha and TH3). Future images with higher reso- 
lution are needed in order to extend the analysis to smaller 
scales. 



6.2 Cross-correlations 

We have introduced a new cross-correlation characteristic 
based on wavelet decomposition, named the wavelet cross- 
correlation, which allows to analyse the correlation as a 
function of scale. The classical cross-correlation coefficient 
is misleading if a bright, extended central region or an ex- 



tended disk exists in the galactic images. In such cases (see 
R, 16, NTH3, CO and DUST in Table 1) the classical cross- 
correlation is dominated by the large-scale structure which 
can be identified as the galaxy's disk with its general radial 
decrease in intensity (see Fig. |2|, bottom right), while the 
(more interesting) correlation on scales of the spiral arms 
can be much worse. In Table 1 this is clearly indicated by 
the wavelet correlation coefficients r„ which are much lower 
than the corresponding values of the classical coefficient rp. 

In classical cross-correlation analysis, the general radial 
decrease of a galaxy's intensity can be removed by dividing 
all values within a certain radial range by its mean (see e.g. 
Hoernes et al. 1998). This is similar to removing the largest 
scale from a map only if the radial decrease is smooth and 
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Figure 11. NGC6946. Wavelet cross-correlations for mid- 
infrared dust emission (DUST) versus total (16), nonthermal 
(NTH3) and thermal (TH3) radio emission. 



independent of azimuthal angle. The wavelet technique is 
recommended as the more reliable method. 

The best correlation, with > 0.75 for any scale 
a, yields the pair DUST-I6 (Fig. The correlation be- 
tween DUST and each of the two components of radio con- 
tinuum, nonthermal and thermal emission, is significantly 
worse. This indicates that the correlation between the total 
intensities is not due to a single process like star-formation 
activity or coupling of magnetic fields to gas clouds, but 
needs a cooperation of physical processes on all scales. A 
detailed analysis will be the subject of a separate paper. 

Similarly, the far-infrared emission of the total dust 
in M31 is highly correlated with the total radio emission, 
while the correlations between the components warm dust 
- thermal radio and between cool dust - nonthermal radio 
are worse (Hoernes et al. 1998; Berkhuijsen, Nieten & Haas 
2000). 

We have shown that the anticorrelation between PIG 
and TH3 (and between PIG and Ha) can be quantified us- 
ing wavelet analysis. This phenomenon, which is due to the 
phase shift between the magnetic arms and the optical arms 
(Beck & Hoernes 1996; Frick et al. 2000), exists only on 
intermediate scales and is clearly visible as a minimum in 
Fig. ^ (lower curve) and Fig. |l^ (right). The position of this 
minimum indicates the scale on which the anticorrelation is 
strongest, i.e. ~ I'.O or l.G kpc at the distance of NGC G94G, 
which is the typical width of a spiral arm. 

The absolute value of the wavelet correlation coefficient 
in these minima is relatively small (0.1-0.25 only). This fact 
can be explained twofold. Firstly, the structures which are 
responsible for this anticorrelation (the arms) occupy only a 
limited part of the image. Secondly, in contrast to the illus- 
trative example given in Fig. H, real arms are more irregular 
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Figure 12. NGC 6946. Wavelet cross-correlation between Ha and 
thermal radio continuum emission (TH3). 



in space and in scales. This is why the position of the local 
extremum in the correlation curve (a) can be informative 
even if the absolute value of r™(a) is relatively small. 

Fig. |l^ shows the correlation between the two tracers of 
ionised gas. The correlation is high on intermediate and large 
scales, but relatively low on small scales. This gives further 
indication for enhanced absorption of the Ha emission on 
small scales. Thus wavelet analysis can be used to study the 
distribution of the absorbing dust. 

Future radio images with better resolution will allow to 
study the correlation on smaller scales. 
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